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Phase-sensitive measurements of the superconducting gap in Fe-based superconductors have 
proven more difficult than originally anticipated. While quasiparticle interference (QPI) mea¬ 
surements based on scanning tunneling spectroscopy are often proposed as defnitive tests of gap 
structure, the analysis typically relies on details of the model employed. Here we point out that 
the temperature dependence of momentum-integrated QPI data can be used to identify gap sign 
changes in a qualitative way, and present an illustration for s± and s++ states in a system with 
typical Fe-pnictide Fermi surface. 


PACS numbers: 74.20.-z, 74.70.Xa, 74.62.En, 74.81.-g 

There is considerable indirect evidence, but currently 
no “smoking gun”, indicating that most, if not all, of 
the Fe-based superconductors condense into a supercon¬ 
ducting state that has Aig symmetry but changes sign 
between electron and hole pockets. This is indeed the 
prediction of spin fluctuation theorjO that follows from 
the simple argument that the pairing interaction, pro¬ 
portional to the magnetic susceptibiltiy, is repulsive and 
peaked at a wave vector Q which nearly nests these pock¬ 
ets. The simplest version of this state, conventionally 
referred to as the s± state, is isotropic on each Fermi 
surface pocket with opposite signs for electron and hole 
sheetfl. The actual gap function realized in Fe-based SC 
is in many cases thought to be highly anisotropic and 
even possess nodes, but may still be considered s± as 
long as the average sign on electron pockets is opposite 
that on hole pockets. 

Arguments against the s± picture have been given as 
well. If orbital fluctuations dominate spin fluctuations, 
a state with equal sign on all pockets, denoted s++, is 
favorecP. In addition, when either hole or electron pock¬ 
ets disappear from the vicinity of the Fermi level, d-wave 
pairing may be enhancec3. Thus, as in the cuprates, the 
determination of the gap symmetry and structure is im¬ 
portant as a clue to the underlying mechanism of super¬ 
conductivity. As discussed in Ref. [H this task is not as 
straightforward as in the cuprates due to the multiplicity 
of Fermi surface sheets. While many tests of s± pairing 
symmetry have been proposed, all seem to rely on spe¬ 
cialized theoretical assumptions or are applicable only to 
certain systems. The phase-sensitive probes that proved 
so decisive in the identification of d-wave symmetry in 
the cuprates are much less useful in the Fe-based super¬ 
conductors, due to the fact that both s± and s++-wave 
symmetries belong to the same Aig irreducible represen¬ 
tation, and the difficulty of fabricating junctions with 
controlled properties. A general qualitative test of an¬ 
other type that could distinguish an s± pair state from 
others would be extremely useful. 


One promising technique in this regard is quasiparti¬ 
cle interference, or Fourier transform scanning tunneling 
microscopy (FT-STM). This probe measures the wave¬ 
lengths of Friedel oscillations caused by disorder present 
in a metallic or superconducting system, which in prin¬ 
ciple contain information on the electronic structure of 
the pure system. These wavelengths appear in the form 
of peaks at particular wavevectors q(a;), which disperse 
with STM bias V = w/e. There is no reliable quanti¬ 
tative theory of quasiparticle interference, in general be¬ 
cause the sources of disorder, exact impurity potentials, 
and k-dependence of the tunneling matrix elements are 
unknown. However, the positions of the q do not depend 
on these effects and are related only to the electronic 
structure, including the superconducting gap function. 

The notion that subsets of these q that connect gaps 
of equal or opposite sign on the Fermi surface can be 
enhanced or not according to the type of disorder was 
introduced by Nunner et al.^^. Pereg-Barnea and Frand^ 
then proposed that a disordered vortex lattice should be¬ 
have like a set of localized order parameter suppression 
scattering centers, and could be used to introduce disor¬ 
der in a controlled fashion with an external field. This 
experiment was performed on the cuprate superconduct¬ 
ing compound, Ca2-xNaa;Cu02Cl2, by Hanaguri et al^, 
who showed that certain quasiparticle interference (QPI) 
q were indeed enhanced or suppressed by the field, in a 
manner apparently consistent with d-wave pairing. 

This experiment was attempted in the Fe-based su¬ 
perconductors on a a Fe(Se,Te) superconducting sample 
near optimal doping by Hanaguri et aP, who identified 
three interband scattering wave vectors q: one associated 
with hole-electron scattering (the smallest momentum), 
and two associated with two different electron-electron 
scatterings. The last two features were enhanced in the 
presence of magnetic field with respect to the first one, 
which led the authors to conclude that the hole and elec¬ 
tron pockets have opposite si gns of the order parameter, 
as expected in the s± modePi^l, Some caveats regard- 
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ing the interpretation, associated with the fact that two 
of the wave vectors coincide with Bragg peaks, were dis¬ 
cussed in Ref. [U and there are other issues that we con¬ 
sider below. More recently, Chi et al.l^i^ measured QPI 
on LiFeAs in zero external field, studying the bias depen¬ 
dence of a large q peak corresponding to interband and 
a small q peak corresponding to intraband scattering, 
arguing that upon moving from large to small bias, the 
fact that one was suppressed while the other enhanced 
could be consistent only with pairing. This is also an 
argument of the qualitative type, but because the actual 
intensity dependence of the peaks as a function of bias 
was in fact nonmonotonic, and because there is currently 
no well-founded theory to support this conclusion, the 
interpretation is uncertain. 

In fact, while QPI measured by FT-STM is in princi¬ 
ple a very powerful qualitative tool, it is nearly useless 
as a quantitative one. This is because, while the poles 
of the response function (q-spot positions) are universal, 
the weights are not, but depend on the type, strength and 
range of the scatterers (not generally known), as well as 
the details of the electronic structure (including the gap), 
and the z-position at which the local density of states 
(LDOS) is calculated. The actual pattern calculated in 
even the most sophisticated theories then generally bears 
little detailed resemblance to the measured STM Fourier 
intensities (for a discussion of some of the causes see 
RefUa), and identifying the symmetry or structure of the 
gap function by a comparison of the intensity patterns is 
similar to determining the age of a model from a Picasso 
cubist painting. 

The situation is further exacerbated by advanced 
mythology that has formed over years about QPI. A part 
of this mythology involves belief that the QPI spectra are 
defined by the geometry of the Fermi surface in a nesting¬ 
like way, another part that in the superconducting state 
these spectra are proportional to coherence factors, and 
that magnetic and nonmagnetic impurities have opposite 
effect on QPI. While this mythology has been partially 
rectified in some papers, it has not become generally ap¬ 
preciated, and some myths still persist. 

A systematic summary of the established facts is badly 
needed, and we shall provide it later in the paper. But, 
the main purpose of this paper is to find a clean quali¬ 
tative way to design a QPI experiment capable of distin¬ 
guishing from other pairings, and in general to identify 
sign-changing gaps in an unconventional superconducting 
system. We argue that the best way is to measure the 
intensity of integrated interband scattering peaks in the 
Fourier transformed density of states, antisymmetrized 
with respect to STM bias, as a function of temperature. 
Then the integrated weights of the set of q correspond¬ 
ing to sign-changing scatterings display a strong enhance¬ 
ments only for this channel; furthermore these qualitative 
distinctions are robust against the strength of the scatter¬ 
ing. We conclude that while the general theory of QPI in 
unconventional superconductors has many pitfalls, this 
particular consequence of pairing sign change is robust 


and can be used to unambiguously identify pairing 
if this type of distinction between intra- and inter-band 
pairing q-peaks is observed. 

I. FORMALISM 

Here we present a model that captures all the qual¬ 
itative features of the general 2-band case, but allows 
analytical evaluation of the density of states. We are 
interested in the density of states in a 2-band supercon¬ 
ductor of isotropic or s± type in the presence of 
disorder of various kinds. Here, following the tradition, 
we use the the words “1-band” and “2-band” for mate¬ 
rials with arbitrary complex normal electronic structure, 
but with the superconducting order parameter A that 
can be approximated by a single value for all states on 
the Fermi level, or by two different values, depending on 
the location in the Brillouin zone. 


A. 1-band problem 

We first remind the reader of the 1-band problem. We 
assume a random distribution of Nj point-like impurities 
at sites with i = 1... Nj. The LDOS can be formally 
decomposed p(r, w) = po{uj) + Sp{r,uj) where po is the 
DOS of the homogeneous superconductor, and 6p is the 
local change in the DOS due to disorder given exactly by 

1 r n 

6p{r,u;) = — Im V G0(r - RQ4 gO(R, - r) 

TT L J 11 

*,1 = 1 

where the w-dependence is suppressed for simplic¬ 
ity, T{u}) is the 2Ni x 2Ni many-impurity T-matrix 
(the factor of two is due to spin), G°(r,a;) is the 
bare advanced electron Green’s function G°(r, w) = 
G’^(k, w) exp(zk • r), " refers to Nambu space, and 
[■ • \ap are Nambu spinor indices. The T-matrix repre¬ 
sents ^the solution to the scattering problem — V + 
V&T for N impurities of identical potential Vi at R^. 
Note that here we have already made an important ap¬ 
proximation, neglecting Umklapp processes (that is, us¬ 
ing one argument in the bare Green’s functions instead 
of two). We will discuss this approximation in the Ap¬ 
pendix A. 

It was shown in Refs. m and M that the poles of the 
many-impurity response (change in LDOS) are identi¬ 
cal to the single-impurity problem, although the weights 
may differ substantially. In addition, in the limit of weak 
impurity potentials, the response for the two problems 
are identical modulo a disorder-dependent factoifi^. For 
clarity, we therefore consider the simpler problem and 
present results for a single impurity, and replace the full 
T-matrix by Tij —>■ iiSij, with T = [1 —V)G°(r = 0)]“^Vi. 

It was shown in Ref. [TH| for simple metals, and gen¬ 
eralized to unconventional superconductors by Ref. 1161 
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that information about the pure electronic system could 
be extracted from STS measurements by examining the 
Fourier transform of dl/dV maps. The Fourier transform 
of the LDOS is p(q, w) = p { y , w ), where 

L X L is a square set of LF' positions at which measure¬ 
ments are made, and q = 27r(m, n )/L are vectors in the 
associated reciprocal lattice. The result for a single band 
is then 


(5p(q,w) = -Im^ [G°(k,a;)t(a;)G°(k-|-q,a;)j (1) 

TT ^^ L J 11 

k 

= ^Tr Im^(To -f r 3 )G'°(k,w)t(a;)G°(k + q.w). 


Here Ti, i = 0,3 are the Pauli matrices spanning Nambu 
space. 


B. 2 band model: q-integrated LDOS peaks 


rately, 

(5p(q~0,a;) (2) 

= ^Tr lm'^{To + T3)Gl{k,uj)i^i,{uj)G°(k + q,uj). 

The form of is known for simple cased^but can be 

a complicated function of the various integrated Green’s 
function components, so we do not specify it here. Full 
expressions are given in the Appendix C. 

Suppose now that we wish to calculate the total weight 
in the small q QPI spot as a function of frequency, defined 
to be 

= ^TrIm ^ (tq- f T 3 )G°(k,a;)fi.^(a;)G°(k + q,w) 

k,q~0,l' 

« ^Tr Im ^ (ro -kr 3 )G°(k,w) 4 ^(w)G°(k-kq,a;) 

k,q,r2 

= ^Trlm {tq + T 3 )Gl{k,uj)i„^{uj)Gl{k',uj), (3) 

k,k'.i/ 


Several works have already establi shed the basic gen¬ 
eralization of § to multiband modelJ^SEBtH] xhe qual¬ 
itatively new aspect is that impurities can scatter quasi¬ 
particles between bands, and if the state is of s± type, 
between bands where the superconducting order parame¬ 
ter Ak has opposite sign. The various types of scattering 
processes can then be classified according to whether they 
connect portions of the Fermi surface with differe nt gap 
sign or not, as was done previously in the cupratej- ^ l ’^ l ^^l 
Without loss of generality for our qualitative purposes, 
we consider a pointlike scatterer with Nambu and band 
space potential V^^{k,k') ~ band in¬ 

dices = ft,, e, where ft and e are just band indices, 
although for typical Fe pnictides they correspond to hole 
and electron pockets, and T 3 and ti correspond, respec¬ 
tively, to nonmagnetic and Andreev scattering, which we 
discuss separately. Weak, purely magnetic scatterers do 
not couple to the spin-averaged local density of states 
measured by a typical STM experiment, so we ignore 
this possibility for the moment (while in Ref. [7] it was in¬ 
correctly suggested that magnetic scattering contributes 
to QPI with a different coherence factor, this mistake was 
corrected in a later paper by the same author^^. Note 
that the generalization to 3 or more bands in the Fe- 
based systems is straightforward and should not change 
our basic conclusions. 

In the Fe-based systems, the Fermi pockets correspond¬ 
ing to bands 1 and 2 are well-separated in momentum 
space and generally have small radius. It is therefore rea¬ 
sonable to expect that isolated spots of scattering inten¬ 
sity will be observed generically at small q, corresponding 
to intraband scatteriiig processes, and at large q, corre¬ 
sponding to interbancpH The intraband term is a simple 
sum of the 1-band expression applied to each band sepa¬ 


where in the second step we extended the sum over the 
small range of q around q = 0 to the full Brillouin zone, 
since the product of two Green’s functions from the same 
band is sharply peaked at small q for small w. In the last 
step, we expressed the double k sum as independent sums 
over the Nambu Green’s functions, which then decouple 
provided that t is momentum independent, as we have 
assumed for the moment. In the simplest approxima¬ 
tion with a flat DOS near the Fermi level, the integrated 
matrix Green’s function is 






WTq -I- AyTl 
\/a;2 - Al 


(4) 


Thus one can use Eq. Q to perform all momentum in¬ 
tegrations in Eq. (§ and obtain closed form expressions 
for the intraband q-integrated LDOS weight 5pintra{<^) 
or the corresponding interband quantity describing scat¬ 
tering between two pockets separated by qo. 


sPinter{u}) = “Trim E (5) 

q~qo 

= ^TrIm E (t-o + T3)G°(k,w)V(a;)G°(k',a;). 

k,k',/j.7^i./ 


Note that we also performed a calculation for the lattice 
(momentum-resolved) based model with two bands, giv¬ 
ing parabolic like electron and hole band dispersions near 
the F and the M points of the Brillouin zone. The results 
are shown in Section IB. Most importantly, the main fea¬ 
tures, especially the T-dependence of the antisymmetric 
correction to the LDOS and its frequency dependence, 
allowing one to distinguish s_|_+ and s± superconducting 
gaps, continue to hold. 
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II. RESULTS: WEAK POTENTIAL 
SCATTERERS 

A. T — 0 frequency dependence 

To complete the solution, the t-matrix for a given im¬ 
purity type must be specified. Here we argue that the 
basic qualitative features of the QPI patterns that are 
sensitive to the sign change of the order parameter (or 
lack thereof) depend only on the Nambu space struc¬ 
ture of the t-matrix, which can be extracted by con¬ 
structing the components of the conductance properly 
symmetrized and antisymmetrized with respect to bias. 
They do not depend on the detailed energy dependence of 
the complex t-matrix, except insofar as impurity bound 
states are created within the gap. Even in this case, the 
question of s± or s++ can be decided by means described 
below. 

We show this by first considering the case of constant t- 
matrix, valid for weak (Born) impurity scattering. While 
in general t has several Nambu components depending on 
the type of scattering, superconducting state and impu¬ 
rity phase shift, it is instructive to focus on one Nambu 
component at a time. For example, if t = t^T^, as e.g. 
for a weak nonmagnetic scatterer, then 


^Pintrai,^') ~ (6) 

2t + T3){ujto + + A^Ti) 

-----—2-A2 -- 


in other words, within this approximation, there is no 
change in the small-q integrated Fourier transform den¬ 
sity of states due to a nonmagnetic intraband scatterer, 
regardless of the relative sign of the two gaps Ah,e- The 
same is not true in general of the interband contribu¬ 
tion Spinterii^) = ^Trlm Sp{q,uj), where qo is the 
wave vector connecting the two Fermi surface pockets. 
In this case we have, again for the T 3 component of the 
t-matrix. 


^Pinter 


(7) 


O 2, T Tr(ro -b T3)(a;To -b A/,Ti)T3(a;To -b AeTi) 
- 27 r tsp/ipeim- 


= -27r tap/jpelm 


- A/, A„ 




which is manifestly nonzero for |Ae| < cu < |Aft|. Fur¬ 
thermore it is easy to show that in the limit when the 
two gap magnitudes become equal, |Ae| —?► |Aft,|, there 
are two distinct cases. For an s± state, sgnA?,Ae < 0, 
SpinterF) = -2FPiP 2 AhS{LJ - |A/i|), while for an s++ 
state, Spinter —t 0. In the more general case with Ah Z 
Ae, the interband response in the s± case remains gener- 
ically much larger than that in the s++ case, with weight 





FIG. 1. Integrated interband density of states 
SPinter/{tspiP2) for constant (weak) T3 t-matrix. Gap 
magnitudes are |Ai|/Tc = 3, |A 2 |/rc = 2, and artificial 
broadening 77 = 10~®. Solid line: s± state; Dashed: s++. All 
other components of Sp are zero for this case. 


concentrated between the two energies ui = |A;i|,|Ae|. 
In Fig. we plot the interband frequency-dependent q- 
integrated LDOS change for a model constant T 3 t-matrix 
to illustrate this difference. Note not only the change in 
sign of the response due to the numerator of Q 
which vanishes at an energy corresponding to the the ge¬ 
ometric mean of the two gaps in this case, but also the 
overall small scale. 

For a realistic scatterer, the Born limit results given 
above are no longer adequate, and the t-matrix aquires 
components in all Nambu channels. It is instructive to 
consider the response of the system to a scatterer with 
constant t-matrix in each of these channels, even if none 
of these cases corresponds to a physical situation with the 
exception of the T 3 (weak potential) scatterer^. Thus 
one can define the experimentally accessible quantities 

pltL (W) = Pi^tr^a (W) ± Pi^tra (-w), 

and use them to make clear qualitative predictions for 
the existence or nonexistence of a strong QPI response 
in the various channels that are independent of the type 
of scatterer, as summarized in Table|^ These features will 
also correspond to peaklike features in the T dependence 
of the integrated QPI intensity, as we discuss below. 

In particular, it is easy to see that the tq component of 
the t-matrix reverses the responses of and com¬ 
pared to the T 3 scatterer. For the intraband part, we get 
immediately Sp^ntraF) = -\^u\) for 

both s± or s_|_+. The interband processes contribute 0 for 
an s± state and — 2 FphPeAhS{uj — Ah) for an s+_|_ state. 
Finally, for a single impurity with constant t-matrix in 
the Ti channel, one easily finds that the intraband sym¬ 
metric densities of states hp^ra singular. The an- 
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tisymmetric interband density of states vanishes 

identically, but the symmetric part can be large. 

The weight is proportional to + ^e), so taking our 

limit |A;i| —>■ |Ae| for qualitative comparison as above, 
we see that the s++ case is singular while the s± case 
vanishes. 

Although a generic nonmagnetic impurity has all 
Nambu components of the t-matrix, which will mix these 
behaviors, they can be isolated to some extent by con¬ 
structing the symmetrized and antisymmetrized densi¬ 
ties of st ates , respectively, as suggested by Maltseva and 
Colemarpl23| This is because the T 3 component of the t- 
matrix generates only odd frequency changes to the den¬ 
sity of states to all orders in perturbation theory, and the 
Tg and Ti component only even ones. The T 2 component 
does not contribute to the change in the Fourier trans¬ 
formed LDOS to all orders in perturbation theory simply 
because of the Trace over Nambu matrices, together with 
the assumption that the gaps are real (i.e. appear with 
Ti). Furthermore, this statement is correct even if the 
particle-hole symmetry is broken (i.e. the normal part of 
the Nambu Green’s functions has additional T 3 compo¬ 
nent). 




intra 

inter 

S++ 


ro,Ti 

To, ri 


5p^~^ 

X 

X 

S± 


ro,ri 

X 


Sp^~^ 

X 

To 


TABLE I. Possibility of singular integrated QPI intensity 
(Fourier transformed density of states) in the symmetric (-I-) 
and antisymmetric (-) channels for s++ and s± superconduc¬ 
tors. Here Ta indicates the presence of a strongly enhanced 
response for an assumed constant t-matrix in the a Nambu 
channel, and the x indicates the absence of one. Magnetic 
impurities have Nambu symmetry and do not modulate 
the total LDOS within the Born approximation (see text). 


B. Thermal average: STM observable 

Our intention is to make clear qualitative predictions 
for observable quantities in STM experiments. Thus far, 
we have shown only the T = 0 results for the artificial 
case of a constant real t-matrix. In any measurement 
at finite temperature, the conductance will be related to 
the change in LDOS Sp convolved with a thermal factor 
weighting the contribution of different electronic states 
to the current. The conventional resulfP^, translated into 
our notation, is 



T/T, 


FIG. 2. Thermally averaged interband q-integrated LDOS 
change within constant density of states and constant t-matrix 
approximation for weak nonmagnetic scatterer with rs (Born 
limit) potential. The external frequency D was taken to be 
2.5 in units where | Ai| and | A 2 I were 3.0 and 2.0, as in Fig. 
Shown is antisymmetrized LDOS {Sp~{Q)) for s++ (dashed) 
and for s± state (solid). Symmetrized components are zero 
for both states. 


in both S 3 __|_ and states in this simple approximation. 
Measurement of the antisymmetric components of the 
interband density of states alone should suffice to quali¬ 
tatively distinguish the two states. Note that the clearest 
results are obtained when the STM bias corresponds to 
an energy between the two gap energies, which can be 
identihed from local tunneling spectra. 

It appears at first glance from Table |T] that measure¬ 
ment of the symmetrized interband q-integrated LDOS 
might also distinguish the two states. However, it is im¬ 
portant to stress here that there is no physical impurity 
in the case of an s-wave superconductor corresponding 
to a situation where the t-matrix is entirely of rg type; 
other components of the f-matrix are mixed at strong 
impurity potentials. This we show below for the realistic 
Coulomb screened impurity potential. In particular, we 
will show that very little difference between s_|_ 3 _ and s± 
will be observed in the symmetrized channel for realistic 
situations. We therefore believe that the symmetrized 
channel should be ignored in the analysis of STM data. 


{Sp^m ^ 

POO 

/ duj6p^{uj) 


( 8 ) 


dw 


-df 

(W -h n) ± - iui — fl) 

OOJ 


In Fig. we plot the nonzero antisymmetrized q- 
integrated density of states for the interband QPI peak 


It is useful to observe that the temperature dependent 
average conductances predicted here bear no relation to 
the standard forms one might expect were quasiparticle 
interference really described by conventional coherence 
factors, as anticipated in Refs. I7I23I This issue is dis¬ 
cussed in some length in Appendix B. 
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III. GENERAL SCATTERING POTENTIALS 


A. Realistic screened Coulomb potentials 


The t-matrix for a single impurity is given by 


t = 


l-C/T3^G(k, 


w 


T 3 U 


(9) 


where the integrated matrix Green’s function given in 
Eq.(4) can be further re-written as ~ 

9 u>,vTo + gA,vT\ for each of the bands. The scattering 
potential matrix V_ can be then separated into intraband 
Uaa — ^bb — Uintra '^0 and the interbaud Uab — 
term. 

We have shown that a constant t-matrix in the Born 
limit leads to a clear prediction of QPI intensities en¬ 
abling one to distinguish from states. For in¬ 
termediate to strong scatterers, however, the t matrix is 
complex and frequency dependent, and includes ti and 
T 2 components in addition to tq and T 3 . One may thus 
be concerned that our conclusions may not be general, 
given that one does not know a priori the strength of 
impurities giving rise to the QPI signal. We therefore 
present results for the full t-matrix of a single impurity of 
arbitrary strength and intra- vs. interband scattering po¬ 
tential, within the flat normal state DOS approximation 
(Eq. 15 . Since for large q a screened Coulomb potential 
falls off like 1 / {q^p -I- q^), and screening lengths 27:/qpF 
are of order the unit cell size, realistic interband scatter¬ 
ing potentials are smaller than intraband potentials. We 
therefore begin in Fig. 0 by hxing Uinter — ^•‘^Uintra for 
various strengths of Umtra- To analyze these results, it 
is useful to identify exactly what terms distinguish s_|__|_ 
and s± states. In the simple pedagogical example with 
which we began, there were no such terms in the intra¬ 
band scattering channel, Eq. , and one term in the 
antisymmetrized interband channel expression for 
Eq. 0 proportional to A/iAg. The full result in the gen¬ 
eral scattering potential case is given in the Appendix C, 
but we find analogously that the intraband LDOS has 
no proportionality to terms sensitive to sign changes. In 
the interband channel there are no A^Ag terms in the 
symmetrized LDOS, whereas the antisymmetrized LDOS 
^Pinter ®till Contains only a single term proportional to 


C/iIm' 


A;,A, 




( 10 ) 


i.e. precisely the same expression as in the simpler ex¬ 
ample of the Born limit. All terms will be multiplied 
by the denominator of the full t-matrix, of course, which 
also contains terms that weakly distinguish and s± 
states but this does not alter our qualitative conclusions. 

Nevertheless, we find that the singular behavior of 
^Pinter ^ ^^e s± case is preserved until the unitary limit 
is reached. In section [III C| below, we exhibit the experi¬ 
mentally observable consequences of this effect. 



w/Tc ^/Tc 


FIG. 3. Integrated density of states components Sp'tntra inter 
for “realistic” screened Coulomb scatterers Uinter = 
0.2Ui„tra, for Uintra = 0.01 (a,b), 0.2 (c,d), and 10 . (e,f). 
Dashed curves correspond to s++ state, solid curves to s±. 
Red and black are even and odd components of the LDOS, 
Sp^ and Sp~ , respectively. Left panels represent intraband, 
right panels interband scattering channels, respectively. 


B. Role of bound states 

In the cases discussed in Fig. [T| no impurity bound 
states are visible in the subgap region below |Ae| = 
2. In general, the formation of impurity bound states 
in a multiband system is more complicated than in a 
1 -band superconductor, and it has been argued that 
such states are indeed nongeneric, requiring fine-tuning 
of the potent ial t o produce a bound state below the 
lower gap edg^^^. In the 2-band model, subgap bound 
states are found only in a very narrow interval around 
a line UinteriUintro) in impurity potential space which 
approaches Uinter = Uintra for strong impuritieffi We 
therefore discuss the case Uinter = Uintra separately here. 

Figures 0a) and 0b) essentially reproduce the results 
of Fig. for a weak scatterer. With increased scatter¬ 
ing strength, however, bound states are formed in the s± 
case, as seen in Fig. 0d). If the bound state is at low 
energy, it steals so much weight from the coherence peak 
LDOS that the hierarchy of intensities represented in Ta¬ 
ble |T] becomes a bit difficult to distinguish. To use QPI 
as a definitive qualitative tool, it may therefore be nec¬ 
essary to consider only systems without impurity bound 
states in the gap, which occurs sometimes, or to mask 
the impurities that give rise to bound states in the spa¬ 
tial window used for the Fourier transform. Of course. 
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FIG. 4. Integrated density of states components &ptntra,inter 
for isotropic scatterers Uintra = Ui„ter for Uintra = 0.01 (a,b), 
0.2 (c,d), and 10. (e,f). Dashed curves correspond to s++ 
state, solid curves to s±. Red and black are even and odd 
components of the LDOS, Sp'^ and Sp~ , respectively. Left 
panels represent intraband, right panels interband scattering 
channels, respectively. 


the observation of bound states in and of itself is strong 
evidence for s± pairing, provided the nonmagnetic char¬ 
acter of the impurities can be reliably assumed. If no 
bound states occur, the hierarchy of LDOS moments of 
Table [I] is clear (compare Fi g. [^ c,d)). 

In the unitary limit. Fig. ERe-f) the bound state has 
already moved through the gap and the basic hierarchy 
is again preserved. Note that for this simple band, this 
limit is essentially achieved already for potentials of order 
Te. 


C. Finite temperatures 

For transparency, we now remove from consideration 
those components of the q-integrated LDOS which are 
not qualitatively affected by a gap sign change, and plot 
in Fig. the thermally averaged, antisymmetrized, in¬ 
terband integrated LDOS for the two states s++ and s±, 
Eq.(8). Here we have used the standard BCS type T- 
dependence for both superconducting gaps. Thermal av¬ 
eraging has the effect of removing many of the sharp 
spectral features at the gap edge, as discussed peda- 
gogically for the case of a constant t-matrix in Fig. 
above. It is clear that the temperature dependent sig¬ 
nal in the s± case is huge and characteristic, whereas in 



FIG. 5. Thermally averaged antisymmetrized interband q- 
integrated LDOS change {5Pinter) nonmagnetic scat¬ 

terers with parameters Uinter = 0-^Uintra, with full t-matrix. 
The external frequency fl/Tc = 2.5, with jAil/Tc and |A 2 |/rc 
taken as 3.0 and 2.0, respectively, (a) refers to the weak 
potential Uintra = 0.01. Shown are antisymmetrized LDOS 
{5p~{Q)) for S++ (dashed) and for s± state (solid), (b) shows 
the same for intermediate strength potential Uintra = 0.2. 
Here we have used the standard BCS type behavior for both 
superconducting gaps. 

the S++ it is small and featureless. Note that the slight 
decrease of {5p~nter) occurs when the width of the peaks 
in the becomes comparable to the thermal broad¬ 

ening of the derivative of the Fermi function. It is also 
clear that the T-dependencies do not resemble the clas¬ 
sic BCS T-dependencies arising from coherence factors 
in, e.g. NMR relaxation and acoustic attenuation cases 
(see the Appendix B). We propose that the measurement 
of {5p~nter) versus temperature is therefore the clearest 
way of identifying a gap sign change using QPI. 

D. Comparison with experiment 

We have shown in the previous section that the tem¬ 
perature dependence of the antisymmetrized integrated 
LDOS is a sensitive measure of the gap sign change in a 
superconductor, although it does not reduce in any limit 
to the BCS-type temperature dependence expected from 
coherence factors. Here we comment on the results of 
Chi et al.ES on LiFeAs, who neither symmetrized their 
integrated conductance maps, nor made use of the sign 
of the LDOS change, but nevertheless obtained what ap¬ 
peared to be a distinct qualitative result. Shown in Fig. 
1^ are the expected \{5pinter,intra){^)\ for both s± and 
S++ states within our model, with the two gap scales 
indicated. While the overall behavior is somewhat com¬ 
plicated, one can see that it is indeed the case, as inferred 
by the authors of Ref. [TTl that the bias dependence for 

an s_|_state immediately below the the upper gap scale 

I A/ll is opposite for the intra- and interband contribu¬ 
tions, while it is the same in the case of an s_|_+ state. 
Note further that had one not known the gap scales ex¬ 
actly, distinguishing between the two states on the basis 
of this type of measurement might have been difficult: the 
relative magnitudes of intra-and interband contributions 
are probably not to be taken seriously, since it is diffi- 
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n/T, 


FIG. 6. Absolute magnitude of thermally averaged unsym¬ 
metrized interband (black) and intraband (red) q-integrated 
LDOS changes \{6pinter,intra}{^)\ for S++ (a) and s± (b) vs. 
Q/Tc for nonmagnetic scatterers with parameters T = 0.2rc 
and Winter = 0.2Uintra, with full t-matrix, and |Ai|/Tc and 
\/S. 2 \/Tc taken as 3.0 and 2.0, respectively, fl = |Ai| and IA 2 I 
are indicated by dashed lines. 


cult to subtract the q weight of the homogeneous system 
from the intraband. The gap scales are presumed to be 
known from STM experiment from direct measurements 
of the coherence peaks in the real space local conduc¬ 
tance oc but they will be shifted somewhat 

from the underlying values by thermal smearing and gap 
anisotropy. The former type of shift is even evident in 
Fig. i In case of more than two distinct gap values 
the results are even more muddled. These caveats are 
among the reasons why we propose that a T-dependent 
measurement at fixed frequency, ideally between the two 
gap scales, should be a more sensitive measure of the gap 
sign change. 

To illustrate what one should expect, we now plot in 
Fig. 0 the antisymmetric component of the interband 
LDOS as a function of bais uj, for different tempera¬ 
ture. This is the quantity that should show the most 
prominent difference between s_|__|_ and s±. In particu¬ 
lar, this component for the case should exhibit a 
sign change at a frequency corresponding to the geomet¬ 
ric mean of the two gap scales, while s± has a finite large 
value there. With further increase of the temperature, 
the main features decrease in amplitude but should re¬ 
main detectable. We expect that this will be the typical 
behavior seen in experiment. 


E. Effect of particle-hole asymmetry: 
momentum-resolved Green’s functions 

To complete our analysis, we also performed a calcula¬ 
tion assuming a lattice based model with two bands, giv¬ 
ing parabolic like electron and hole band dispersions near 
the r and the M points of the Brillouin zone. In partic¬ 
ular, we consider the simplest band topology of the iron- 
based superconductors with hole band, ^^.(k) = — 

centered near the T-point of the Brillouin zone and the 

electron band, £e(k) = -Mei centered near the 

M point of the Brillouin Zone (Q= (7r,7r)). Here, we 
set ^ ^ = 34 and pe = -10.6, ph = 9.4 (all in 



FIG. 7. Thermally averaged LDOS changes for antisymmet¬ 
ric, interband channel, \{5pinter,intro) (yt)\ vs. fl/Tc. Gurves 
shown are for nonmagnetic scatterers with parameters T/Te = 
0.1 (a), 0.3 (b), 0.6 (c), 0.9 (d) and Winter = Q.2Wintra, with 
full t-matrix, and \/S.\\/Te and |A 2 |/rc taken as 3.0 and 2.0, 
respectively. D = |Ai| and IA 2 I are indicated by dashed lines. 
Solid curve: s±. Dashed curve: s++. 
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FIG. 8. Integrated density of states components 
for isotropic scatterers Winter = 0.2Wintra for Wintra = 0.01 
(a,b), 0.2 (c,d), and 10 (e,f) (all in units of energy as in for 
the momentum independent model), computed for momen¬ 
tum resolved Green’s functions, Eq.(25). As in Fig.3, the 
dashed curves correspond to s++ state, solid curves to s±. 
Red and black are even and odd components of the LDOS, 
and 5p~ , respectively. Left panels represent intraband, 
right panels interband scattering channels, respectively. 
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FIG. 9. Thermally averaged antisymmetrized interband q- 
integrated LDOS change {5{0) for nonmagnetic scat- 
terers with parameters Uinter = O-^Uintra, with full t-matrix 
and momentum resolved Green’s function, Eq(25). The ex¬ 
ternal frequency fl/Tc = 2.5, with |Ai|/Tc and |A 2 |/Tc taken 
as 3.0 and 2.0, respectively, (a) refers to the intermediate 
potential Uintra = 0.2. Shown are antisymmetrized LDOS 
{5p~{Q)} for S++ (dashed) and for s± state (solid), (b) shows 
the same for the larger strength potential Uintra = 10. As in 
Fig.5, here we have used the standard BCS type behavior for 
both superconducting gaps. 



n/T, n/T^ 


FIG. 10. Thermally averaged LDOS changes for antisym¬ 
metric, interband channel, \{5pinter,intra){yi')\ vs. fl/Tc- 
Gurves shown are for nonmagnetic scatterers with parame¬ 
ters T/Tc = 0.1 (a), 0.3 (b), 0.6 (c), 0.9 (d) and Uinter = 
0.2Uintra (Uintra = 0.2), with full t-matrix, momentum- 
resolved Green’s functions, Eq.(25), and |Ai|/rc and |A 2 |/rc 
taken as 3.0 and 2.0, respectively, fl = |Ai| and IA 2 I are 
indicated by dashed lines. Solid curve: s±. Dashed curve: 

S++. 


the same units of energy as in the main text). The most 
important part introduced by the lattice based models is 
the electron hole asymmetry in the normal part of the 
Nambu Green’s functions. 


G°(k,zw„) ~- 


iuJnTo + £;.(k)r3 -I- A^ti 
W^+A2-b£2(k) 


( 11 ) 


Substituting the Green’s function for the electron and 
hole bands, Eq.(ll) into Eqs.(3),(5) and (9) we com¬ 


puted the interband and inraband corrections to the 
LDOS similar to the main text. In Eigs|8][T0| we show 
the corresponding results for the integrated density 
of states components ^pfntra inters thermally averaged 
antisymmetrized interband q-integrated LDOS change 
{Sp~nter) thermally averaged LDOS changes for 

antisymmetric, interband channel, | {Spinter,intra )(D)| vs. 
il/Tc, respectively. Despite some differences, the main 
features that allow one to distinguish s++ and s± super¬ 
conducting gaps continue to hold also if one uses more 
realistic Green’s functions. We believe this happens be¬ 
cause the electron and hole Fermi surfaces and corre¬ 
sponding gaps are well separated in q space. This allows 
one to disentangle intra- and inter-band scattering pro¬ 
cesses in the normal and superconducting states clearly, 
as found in most of the ferropnictides. In particular, we 
observe that the actual T/Tc dependence for the s±-wave 
symmetry shows a non-monotonic dependence which is 
reflected in the peak-like structure at small T/Tc val¬ 
ues and corresponding downward behavior at low T/Tc 
values. This is absent for the s++-wave symmetry. In 
addition, we see from Fig. that the antisymmetrized 
STM signal for |Ai| < 17 < IA 2 I has one sign for s± and 
changes sign for s++. Therefore we believe that these 
two key features can be observed in the experiments. 


F. Remarks on magnetic and ri (“vortex”) 
scatterers 

A weak “magnetic impurity” represented by an iso¬ 
lated classical spin that couples via exchange to con¬ 
duction electron spin density may be written in Nambu 
space as VspinToO's, where is the Pauli matrix in spin 
space. Within the Born approximation, the change to the 
up-spin LDOS will cancel that of the down-spin LDOS. 
Higher order magnetic scatterings will produce effects, 
but since in general a chemical substituent with a mag¬ 
netic moment will have a nonmagnetic scattering poten¬ 
tial much larger than its magnetic one, we ignore this ef¬ 
fect here. Including transverse spin couplings (or deep d- 
or /-levels within an Anderson model approach) will re¬ 
sult in Kondo physics which obviously produces interest¬ 
ing effects on the density of states, including Kondo res¬ 
onances near the Fermi level, with concomitant influence 
on QPI; these have been recently discussed elsewher^^. 
Similar effects on the QPI would be expected when the 
Yu-Shiba bound state is induced by the magnetic impu¬ 
rity. In order to draw qualitative conclusions regarding 
gap symmetry using the methods described here, samples 
or regions of samples displaying Kondo and/or Yu-Shiba 
bound state resonances should be excluded. 

Scattering in the ti, or Andreev channel has been dis¬ 
cussed in several contexts in the field of unconventional 
superconductivity. Ghemical impurities suppress the or¬ 
der parameter in their vicinity, creating an effective off- 
diagonal local pote ntial w hich contributes to the scatter¬ 
ing of quasiparticled291MI Normally these effects are ig- 
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nored, e.g. in standard t matrix calculations, or treated 
as weak, but under some circumstances they can become 
important. If an impurity has the effect of enhancing the 
pairing locally, as occurs in some modelPi^t^, the ti po¬ 
tential component of an impurity can be significant and 
even control the behavior of the conductance spectrum 
and mafP^. 

The order parameter is also suppressed near vortex 
cores, and Pereg-Barnea and Franz suggested that this 
fact could be used to provide a method of controlling dis¬ 
order and distinguishing gap symmetries in situ, provided 
the vortex lattice were sufficiently disordered^. Here we 
do not discuss ti chemical impurities in detail, as we are 
focussed primarily on qualitative aspects of QPI, but we 
discuss the oft-repeated statement that the effect of the 
disordered vortex lattice, represented by a random, tun¬ 
able set of Ti scatterers can distinguish sign- preser ving 
and sign-reversing QPI peaks. It is belie vecPI^Ei^ that 
the peaks whose weight is enhanced in a field correspond 
to sign-preserving peaks, while those whose weight is sup¬ 
pressed by a field are sign-reversing. 

In fact a clear statement to this effect is difficult to 
make. As we showed in Sec. |II A[ ri impurities indeed 
enhance the QPI signal in the intraband channel, which 
represents sign-preserving scattering for both supercon¬ 
ducting states considered. In addition, Table|T]also shows 
that the sign-preserving interband scatterings in the s++ 
case give rise to an enhancement. There is no indica¬ 
tion of a suppression with field in the sign-reversing (s± 
interband) case, however. This apparent discrepancy 
was noted already by Pereg-Barnea and Frand^, who sug¬ 
gested that the disordered vortex lattice led to a ran¬ 
dom phase potential experienced by quasiparticles, which 
might give rise to an overall suppression of the QPI sig¬ 
nal. Such a “background” suppression could then be 
overcompensated by the singular enhancements of the 
sign-preserving scattering wave vectors. There is, how¬ 
ever, no calculation to support this assertion, so state¬ 
ments about determining gap symmetry from QPI peaks 
in unconventional superconductors from their field be¬ 
havior should be treated with caution (see also the dis¬ 
cussion in the Appendix). 


G. Application to other states; nodeless d-wave 

While for pedagogical reasons we have restricted our¬ 
selves to two bands, isotropic gaps, and either or s± 
states, the concepts we have discussed are clearly appli¬ 
cable to more general situations. The interband entries 
of Table |l] labelled s++ apply generally to gap sign pre¬ 
serving transitions between bands, and those labelled s± 
apply to sign-changing ones. The obvious example under 
discussion in the Fe-based superconductivity field is the 
putative d wave state in systems with no hole pockets but 
four electron pockets at the (±7r, 0) and (0, ±7r), int itially 
proposed for the alkali-intercalated FeSe material d^^ l ^^l 
In addition to Bragg scattering and small q intraband 


scattering, interband scattering should give rise to rings 
of scattering intensity around q = (0, ±27r) and (±27r, 0), 
as well as around (±27r, ±27r). Predictions for intraband 
scattering weights will be identical to those listed for s± 
in Table m 


IV. CONCLUSIONS 

In this paper we have argued that the task of identi¬ 
fying order parameter symmetries in unconventional su¬ 
perconductors via QPI measurements is unlikely to be 
successful if it relies, as in the past, on comparisons of 
theoretical conductance maps with experiment. This is 
because there are too many unknown parameters, par¬ 
ticularly in multiband systems, to allow for a quantita¬ 
tive theory of QPI. It is possible that this situation can 
be improved to some extent by ab initio based calcula¬ 
tion of the density of states away from the metal surface, 
including an isolated impurity. This would be however 
time consuming and would still not eliminate all quanti¬ 
tative assumptions. On the other hand, in systems like 
the FeSC, where one can clearly identify QPI peaks re¬ 
lated to intra- and interband scattering, the temperature 
dependence of the integrated weights of these peaks can 
provide a robust qualitative means of detecting order pa¬ 
rameter sign changes that can be used to determine its 
structure. 

We have focussed our attention on a model of a 2- 
band system with two distinct gaps A, and shown that 
the most sensitive way to distinguish scattering processes 
connecting gaps of same or different sign is to operate 
with STM bias in the energy region between the two gap 
scales, identihcation of which is a relatively simple exper¬ 
imental task. We have then shown that the T-dependent 
response of the symmetrized and anitsymmetrized combi¬ 
nations of the conductance for both intra- and interband 
scattering provide a characteristic signature of a gap sign 
change or lack thereof. These temperature dependencies 
do not, as suggested by previous works, correspond to 
thermal averages of simple BCS coherence factors, but 
are somewhat more complicated. In particular, we find 
that the effect is strongest at low temperatures, and not 
near Tc, in contrast to the expectation assuming coher¬ 
ence factors. Although we have focussed on the question 
of distinguishing and s++ states in the FeSC, it is clear 
that similar arguments can be made for sign-changing 
gaps in other contexts, for instance a putative d-wave 
state in FeSC materials with electron pockets only. 

In the past, most QPI experiments have focussed on 
the power spectrum of the LDOS, in other words the ab¬ 
solute magnitude of the density of states Fourier trans¬ 
form |dp(q)| or related ratios of this quantity, so-called Z 
or R maps. We have shown here that measurement of the 
signed symmetrized and antisymmetrized QPI maps are 
crucial to extract symmetry information in the absence 
of detailed knowledge of the impurity potentials, which 
is usually the case. This effect persists up to Tc. 
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The existence of order parameter bound states is one 
aspect which must be treated with care in such a mea¬ 
surement, as we have shown that they tend to steal 
weight from the spectral region where the characteris¬ 
tic distinctions are most visible. Of course if one has 
a clear indication of a bound state induced by a non¬ 
magnetic impurity, it is already a strong indication of a 
sign-changing order parameter. Nevertheless additional 
complementary evidence can be obtained by performing 
the analysis suggested here while masking the regions 
containing the bound states before Fourier transforming. 

Finally, we have discussed the commonly used method 
of distinguishing gap symmetries by observing the mag¬ 
netic field dependence of QPI peaks and identifying sign¬ 
preserving scattering wave vectors as those correspond¬ 
ing to peaks that increase with field, and sign-changing 
ones with those that decrease with field. This analysis, 
while appealing and possibly correct in some cases, is 
based on a questionable analogy of vortices as pointlike 
Andreev scatterers that may fail for several reasons, in¬ 


cluding if the vortex lattice is too ordered or coherence 
lengths too large. In addition, while we can understand 
within the work presented here why some QPI peaks can 
be enhanced by the magnetic field, there is no firm theo¬ 
retical ground for interpretation of those peaks that are 
suppressed. 

Acknowledgements. We thank P. Coleman, D.J. 
Scalapino and I. Vekhter for useful discussions. P.J.H. 
was supported by NSF-DMR-I005625, and I.I.M. by the 
U.S. Office of Naval Research through the Naval Research 
Laboratory’s Basic Research Program. The work of DA 
and IE was supported by the Focus Program 1458 Eisen- 
Pniktide of the DFG, and by the German Academic 
Exchange Service (DAAD PPP USA no. 57051534). 
P.J.H. and I.I.M. would like to thank for hospitality R. 
Valenti and the Goethe University of Frankfurt, where 
this project was started, as well as the Kavli Institute for 
Theoretical Physics, where several key discussions took 
place. I.E. acknowledges the support allocated to Kazan 
Eederal University for the project part of the state as¬ 
signment in the sphere of scientific activities. 


^ P. J. Hirschfeld, M. M. Korshunov, and I. I. Mazin, Re¬ 
ports on Progress in Physics 74, 124508 (2011) 

^ 1. 1. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008) 

® H. Kontani and S. Onari, Phys. Rev. Lett. 104, 157001 

( 2010 ) 

S. Maiti, M. M. Korshunov, T. A. Maier, P. J. Hirschfeld, 
and A. V. Chubukov, Phys. Rev. Lett. 107, 147002 (2011) 

® T. S. Nunner, W. Chen, B. M. Andersen, A. Melikyan, 
and P. J. Hirschfeld, Phys. Rev. B 73, 104511 (2006), 

® T. Pereg-Barnea and M. Franz, Phys. Rev. B 78, 020509 
(2008) 

' T. Hanaguri, Y. Kohsaka, M. Ono, M. Maltseva, 

P. Coleman, 1. Yamada, M. Azuma, M. Takano, 

K. Ohishi, and H. Takagi, Science 323, 923 (2009) 
http://www.sciencemag.org/content/323/5916/923. full.pdf 

® T. Hanaguri, S. Niitaka, K. Kuroki, 
and H. Takagi, Science 328, 474 (2010) 

http://www.sciencemag.org/content/328/5977/474. full, pdf, 

® F. Wang, H. Zhai, and D. Lee, Phys. Rev. B 81, 184512 

( 2010 ). 

A. Akbari, J. Knolle, 1. Eremin, and R. Moessner, Phys. 
Rev. B 82, 224506 (2010) 

S. Chi, S. Johnston, G. Levy, S. Grothe, R. Szedlak, 

B. Ludbrook, R. Liang, P. Dosanjh, S. A. Burke, A. Dam- 
ascelli, D. A. Bonn, W. N. Hardy, and Y. Pennec, Phys. 
Rev. B 89, 104522 (2014) 

A. Kreisel, P. Choubey, T. Berlijn, W. Ku, B. M. Andersen, 
and P. J. Hirschfeld, Phys. Rev. Lett. 114, 217002 (2015) 

L. Capriotti, D. J. Scalapino, and R. D. Sedgewick, Phys. 
Rev. B 68, 014508 (2003) 

L. Zhu, W. A. Atkinson, and P. J. Hirschfeld, Phys. Rev. 

B 69, 060503 (2004) 

P. T. Sprunger, L. Petersen, E. W. Plummer, E. Lgs- 
gaard, and F. Besenbacher, Science 275, 1764 (1997) 
http://www.sciencemag.org/content/275/5307/1764. full, pdf 
J. E. Hoffman, K. McElroy, D. H. Lee, K. M. Lang, 


H. Eisaki, S. Uchida, and J. C. Davis, Science New Se¬ 
ries, 297, pp. 1148 (2002) 

S. Sykora and P. Coleman, Phys. Rev. B 84, 054501 (2011) 

T. Hanke, S. Sykora, R. Schlegel, D. Baumann, L. Har- 
nagea, S. Wurmehl, M. Daghofer, B. Buchner, J. van den 
Brink, and C. Hess, Phys. Rev. Lett. 108, 127001 (2012) 
F. Wang, H. Zhai, and D. H. Lee, EPL 85, 37005 (2009). 
E. Plamadeala, T. Pereg-Barnea, and G. Refael, Phys. 
Rev. B 81, 134513 (2010). 

Y. Y. Zhang and et ah, Phys. Rev. B 80, 094528 (2010). 

Q. -H. Wang and D.-H. Lee, Phys. Rev. B 67, 020511 
(2003) 

M. Maltseva and P. Coleman, Phys. Rev. B 80, 144514 
(2009) 

D. V. Efremov, M. M. Korshunov, O. V. Dolgov, A. A. 
Golubov, and P. J. Hirschfeld, Phys. Rev. B 84, 180512 
( 2011 ) 

In a d-wave superconductor, the large impurity potential 
limit yields t oc tq, but in an s-wave system the nonvan¬ 
ishing integrated anomalous Green’s functions generate ri 
and T 2 terms as well. 

J. E. Hoffman, Reports on Progress in Physics 74, 124513 

( 2011 ) 

R. Beaird, I. Vekhter, and J.-X. Zhu, Phys. Rev. B 86, 
140507 (2012) 

M. A. Derry, P.G. and D. Logan, arXiv: 1503.04712 
(2015). 

M. H. Hettler and P. J. Hirschfeld, Phys. Rev. B 59, 9606 
(1999) 

A. Shnirman, i. d. 1. Adagideli, P. M. Goldbart, and 
A. Yazdani, Phys. Rev. B 60, 7517 (1999) 

I. N. Khlyustikov and A. 1. Buzdin, Adv. Phys. 36, 271 
(1987). 

T. S. Nunner, B. M. Andersen, A. Melikyan, and P. J. 
Hirschfeld, Phys. Rev. Lett. 95, 177003 (2005) 

M. M. Maska, i. d. Z. Sledz, K. Czajka, and M. Mierze- 






12 


jewski, Phys. Rev. Lett. 99, 147006 (2007) 

K. Foyevtsova, R. Valenti, and P. J. Hirschfeld, Phys. Rev. 
B 79, 144424 (2009) 

K. Foyevtsova, H. C. Kandpal, H. O. Jeschke, S. Graser, 
H.-P. Cheng, R. Valenti, and P. J. Hirschfeld, Phys. Rev. 
B 82, 054514 (2010) 

M. Yang, F.and Gao, Z.-Y. Lu, T. Xiang, and D.-H. Lee, 
Europhys. Lett. 93, 57003 (2011). 


T. Maier, S. Graser, P. Hirschfeld, and D. Scalapino, Phys. 
Rev. B 83, 100515 (2011). 

Z. P. Yin, K. Haule, and G. Kotliar, Nature Phys. 10, 845 
(2014). 

F. Ahn, 1. Eremin, J. Knolle, V. Zabolotnyy, S. Borisenko, 
B. Buchner, and A. Chubukov, Phys. Rev. B 89, 144513 
(2014). 

M. Tinkham, Introduction to Superconductivity: Second 
Edition , Dover Books in Physics (2014). 


V. APPENDIX 


A. Remarks on the Bragg peaks and the 1-Fe vs. 2-Fe Brillouin zone 

The electronic structure of the iron-based superconductors can be described exactly in the Brillouin zone corre¬ 
sponding to a 2-Fe unit cell, or approximately in the twice larger one, corresponding to a 1-Fe celffl In the former 
case, the intraband scattering create a QPI spot near the zone center (q « 0), and an intraband one near zone cor¬ 
ners (q «7r/a, tt/o), where a is the lattice parameter of the 2-Fe cell, a = dV^, and d is the Fe-Fe bond length. In 
the experiment, the QPI signal is invariably mixed with the Bragg peaks, resulting from electrons scattering off the 
regular crystal lattice. It is important to understand that, just as in X-ray scattering, while the density of states p(q) 
is peaked at each reciprocal lattice G, the intensity of each q -I- G component depends on G, and generally decays 
with |G|, and the same is true for the QPI spectra. In other word, when we say that the intraband spot is located 
around q = 0, it is implied that there are also spots near all q == G, albeit with a reduced intensity. 

Now, let us use the 1-Fe cell, and a twice larger Brillouin zone. In this setting, there are spots near qi« 0, near 
q2«7r/d, 0, and two overlapping spots near q3_4ft!7r/d, n/d. Three of them, the first one and the last two, are located at 
2-Fe reciprocal lattice vectors. Of course, this means that as long as we include scattering off the pnictogen or chalcogen 
sublattice, each of these peaks will generate shadow peaks at all other reciprocal lattice vectors. Experimentally, 
however, they will be clearly distinguishable: the qi peak will be most intensive at G = 0, while qa^ will actually be 
stronger at G = 0,±7r/a and ±7r/a, 0 (in the reduced zone). The spots near qa will not be affected by the downfolding 
procedure, as for this particular vector the sublattice scattering will not generate any shadows. 

This may seem to be of academic importance, but it may have considerable practical ramifications. Indeed, while for 
Fe-based superconductors, as mentioned, there is no problem separating the hole-hole and electron-electron scattering 
from the hole-electron one, the scattering between two inequivalent electron pockets in the 2-Fe zone is seemingly 
indistinguishable from the intraband scattering, and, as the latter, overlaps with the Bragg peaks. But, as discussed 
above, the intraband scattering will be stronger at that half of the reciprocal lattice vectors that coincide with those 
reciprocal lattice vectors of the 1-Fe cell, while the interband electron-electron scattering will be stronger at the other 
half of the G-vectors. This, of course, makes the signatures of the order parameter signs discussed in the main text 
weaker, but does not destroy them. Mo reover , several papers recently suggested a possibility of the order parameter 
sign change between different hole band^SHMl particular, between the two bands formed by the xz and yz orbitals 
and the one with predominantly xy character. In the 2-Fe zone they both occur at the zone center and thus seem 
indistinguishable. However, a closer look reveals that the first two bands generate a spot that is located at the zone 
center in both 1-Fe and 2-Fe zones, while the scattering between xzjyz and xy pockets, in the 1-Fe zone occurs near 
the zone corner. While both processes will create “shadows” at all 2-Fe reciprocal lattice vectors, the former will be 
stronger at one half of the vectors, and the other at the other half, again allowing to use QPI to assess the above 
hypothesis. 

Finally, the very fact that some of the nontrivial QPI spots are overlapping with Bragg peaks is not a disadvantage, 
but just the opposite. As pointed out by Hanaguri et aL“ the Bragg peaks are much sharper then the QPI maxima, 
and have distinctly different profiles, which allows them to be separated. Note that the mechanism that is supposed 
to generate a QPI dependence on the external magnetic field should not operate on the Bragg peaks, whether in the 
normal or superconducting state. Thus, if not only the QPI, but also the Bragg peaks demonstrate a strong field 
dependence (as was in fact the case in Ref.^^, this strongly suggests that the field dependence of the QPI spectra 
is also not directly related to scattering off Abrikosov vortices. Given the caveats described in the main text, this 
capability appears rather useful. 
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B. Coherence factors 


We remarked in the introduction that there is an expectation shared by many in the STM community that the 
temperature and bias dependence of QPI signals will follow simple BCS coherence factors. This is based in part on 
Refs. I7I23I which expressed the change in the Fourier transformed density of states as an integral over such factors, 
written explicitly, for instance, for the antisymmetrized response and Born scattering, as (u^Uk' and similar 

expressions for other types of scattering. However, this assertion is not correct and the prefactors in the expressions 
for the QPI intensity cannot be cast into such form. Let us illustrate that now for the Born scattering. In this case, 


Sp(q,uj) = ^Trlmy r3Gk(a;)r3G'k'(w)< 5 (k - k'-q) 
f + ^k^k+q - AkAk' ^ ^ 


( 12 ) 

(13) 


where ^ are the one-electron energies. Concentrating on the fraction under the integral, we see that it is 

— AkAk 
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(u;^-E^) 


S(uj — Ek)(E^ -I- 44' — AkAk') 


4Ek(E^-E^,) 

S(uj — Ek’)(E^, + 46k' — AkAk') 
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1 ,S(u;-Ek)(EE 


■ 6k6k' — AkAk' 




4Ek 

S(uj — Ek>)(Ek^ + 6k6k' — AkAi, 


4Ek^ 


(14) 


Note that arguments of the <5—functions are different therefore in order to combine the two terms we need to rename 
the k variables in the second term, after which, assuming inversion symmetry in A and 6, we get: 


Sp(q,uj) = 


S{u! — Ek){E^ -I- 6k6k' — AkAk') 


Ik ■- 


8Ek{Ek^ - El,) 


S{uj — EkXjEl, -I- 6k6k' — AkAk') 

4EkXEl, - El) 


5(k — k'—q) 


{El + 6k6k' - AkAk')4c'^ - Ek) <5(k - k'-q) 
{El-El,)u: 2 


(15) 


The first factor is what was taken to be a coherence factor in Refs. [7|and[23l However, the true coherence factor in 


question is different, namely {El + 6k6k' ~ AkAk')/(i?k£'k')j that instead of the expression above one would have a 


very different formula, namely 




f (^^ + 6k6k' - AkAk')<I(a; - Ek)S{u; - Ek') (5(k - k'-q) 

k 2 

(a;2 + 6k6k' - AkAk')<5(w - Ek) ^(k - k'-q) 


ui{u;^-El,) 


(16) 


One may think that this difference will disappear after an integration over q, and the integrated expression will be 
similar to classic BCS predictions of thermally averaged transition probabilities, e.g. in NMR spin relaxation and/or 
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ultrasonic attenuatioiP^. For example, consider the spin-lattice relaxation rate 

^ „ ii„ v 

TiT ljq—J-O Cu’q 

~ lim —Im Tr (dG'°(k,w„)Q;G°(k',a;„-I-Wo)) , 

cJo—^■O CJq ^ / iuJo^uJo-\-iO+ 

\c \c' II 7V 


(17) 


where x is the spin susceptibility matrix in band space, a is the electronic spin operator in Nambu space, wq is the 
Larmor frequency, and we have neglected any momentum dependence of hyperfine matrix elements in order to arrive 
at a simple expression. In fact, Eq. 0 bears a certain resemblance to the structure of Eq. Q , in particular because 
each Green’s function is local and therefore integrated over momentum independently. Performing these sums and 
the analytical continuation, one obtains the usual result, with slight generalization to multiband systems: 


1 





= intraband terms -I- 2 



2 






Im 




f'OQ 

= intraband terms + 2 duj 


> max(A/i,Ae) 


df\ up- + AhAe 

duo) - Al^uj^ - Al 


(18) 


Now one can see that an s_|_+ superconductor will have an interband contribution to NMR relaxation that obeys 
exactly a BCS type II coherence factor T-dependence (in the limit A/j —>■ Ae), while for s± the corresponding result 
is type I. 

Compare now to our expression for the interband contribution to the Fourier transform density of states in the case 
of a weak potential scatterer, Eq. Q, together with the thermal average. 


{^Pinter{uo')') CX 



df '\ - AftAe 

duo) A2y/w2- A2 

^_a/\ - A,Ae 

V duo) ^ A^' — uo"^ — A2 


(19) 


We see that while the functional form of the fraction under the integral is indeed the same, the cost of the q— integration 
is that the frequency integral is now taken over a totally different, in fact not overlapping, limits, which results in a 
completely different T-dependence, as shown above. 


C. Single impurity t-matrix for 2-band system 

The t-matrix for a single impurity is given by 


t = 


1 - ]Lt'a X] S(k, w) 

k 


r-iU_ 


( 20 ) 


where the integrated matrix Green’s function given in Eq.Q can be further re-written as (X^k S:(kj ^))y = goj,iy'’'o + 
5 a,i/T for each of the band. The scattering matrix G can be then separated into intraband {iQ^v = Uintra ^^nd 
interband {iQ^v = Umter term. 
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Inverting the t matrix one finds the common denominator for both intra and interband terms of the t matrix 

^ i,9uj,h 9A,h) ^intra i,9uj,e 9A,e) ^intra ^ 9A,h9A^e) ^inter 

+ {9t,h ~ 9 a, h) ( 5^,2 “ 5 a , 2 ) {Winter ~ Sintra) ^ + 1 - 
Now for the intraband scattering within the band 1 one obtains in Nambu components 

ih = t^To + tin + tin 


tfi — [5<^,?it^mtra + 9ui,eUj,nter 9u>,h {9u>,e 9A,e) {^inte 

[ 9A,hUintra 9A,eUinter 9A,h {^9ixj,e 5A,e) 

~ [Uintra (l ~ (5cj,e ~ 9A,e) iPintra ~ ^inter))] / 


—U^ ] 

inter intra) 


inter ^ intra 


/D, 

■a) /D, 


The expressions for ie are obtained by exchanging the band indices h and e. Similarly, the interband t-matrix can be 
written as 

ieh = tlhTo + tlhTi + tlhT2 + tl^n (24) 


^eh [i9u^,e + UijitraUinter\ /^ 

^eh ~ [ i9A,e H“ 5A,/i) UintraUinter\ /^ 

^eh ~ [~* {9A,e9ai,h ~ 9uj,e9A,h) Uinter (t^inter ~ ^mtro)] 

^eh ~ pinter {9i^,e9oj,h ~ 9A,e9A,h) Pintra ~ ^inter^ 4” l] / 


and the corresponding expression for the- Note that the interband scattering terms generate a T 2 contribution for the 
interband t-matrix. 



